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We present a segregated strain model that can describe the thickness-dependent dielectric proper- 
ties of epitaxial ferroelectric films. Using a phenomenological Landau approach, we present results 
for two specific materials, making comparison with experiment and with first-principles calculations 
whenever possible. We also suggest a "smoking gun" benchtop probe to test our model. 



I. INTRODUCTION 

Increasing demands for high-density data storage with 
ultra-fast accessibility present tremendous challenges. In 
parallel to the characterization of new materials, impor- 
tant size-dependent effects must be understood to op- 
timize design. This is particularly true for ferroelec- 
tric memories, whose nonvolatility and low power con- 
sumption make them well-suited for portable applica- 
tions;^'^ their dielectric properties are strongly depen- 
dent on electromechanical boundary conditions due to 
the long-range nature of their underlying electrostatic in- 
teractions. The sensitivity of ferroelectricity to homoge- 
neous strain in bulk perovskite oxides is well-known.- In 
thin films, the effects of homogeneous misfit strain have 
been identified)^ studied and controlled to the point that 
particular systems have been strain-engineered to have 
spontaneous polarizations significantly larger than those 
in the bulki^i^ Despite these impressive achievements, 
several authors^i^iSiiiS have suggested that homogenous 
epitaxial strain cannot qualitatively account for all the 
observed effects in ferroelectric films. In particular, the 
thickness-dependence of their dielectric properties has 
been attributed to defect-induced strain gradients.^-'- In 
this paper we explore whether these observed size-effects 
are also consistent with an alternative scenario where 
the majority of each ferroelectric film is homogeneously 
strained. Our phenomenological study indicates that 
thickness-dependent dielectric measurements are insuf- 
ficient to determine the presence/ absence of underlying 
inhomogeneous strain and we suggest further benchtop 
experiments that will resolve this issue. 

In well-screened coherent epitaxial ferroelectric thin 
films, uniform polarization is energetically favored. Lat- 
tice mismatch between the film and the substrate is a 
key source of macroscopic strain in these systemsi^iii 
For very thin films, the energy cost for homogeneously 
straining the film to match the substrate is less than that 
associated with the introduction of stress-relieving de- 
fect structures at the interface. However in thicker films, 
such misfit dislocations form and produce strain gradi- 
ents^ i^^'^^i^^'^'* ; threading dislocations and point defects 
are additional sources of inhomogeneous strain. In pla- 
nar ferroelectric films, inhomogeneous strain can affect 
the ferroelectric transition through both smearing and 
through its coupling to the polarization. Such flexoelec- 
tric coupling of strain gradients to the polarization has 
been the topic of much recent interesli^ i^^'^^ particularly 



as it has been suggested that flexoelectric effects are en- 
hanced by large dielectric coefficients Recently it has 
been argued that such strain/stress gradients are crucial 
for the modelling of thickness-dependent dielectric prop- 
erties of ferroelectric films jii^ Here we propose an alter- 
native model: that the misfit dislocations reside within 
a thin buffer layer next to the interface; the majority of 
the film, which is relatively defect-free, is then homoge- 
neously strained. "'^^ In our phenomenological treatment, 
we also include a bulk anisotropyiS in the form of an 
effective field, possibly due to asymmetry of the elec- 
trodes and/or to the thin buffer layer. We model the 
thickness-dependent dielectric properties in two different 
types of ferroelectric films, and compare our results with 
experiment and with first-principles calculations when- 
ever possible. Finally we discuss a benchtop "smoking 
gun" probe to distinguish our segregated strain scenario 
with that of inhomogeneous strain in ferroelectric thin 
films. 

The structure of this paper is as follows. In Section 
II we review the experiments that motivate this study 
and their implications for any descriptive model. Details 
of our phenomenological Landau approach are presented 
in Section III, with specific discussion of the appropri- 
ate boundary conditions and depolarization effects. In 
Section IV we present our main results for films of two 
specific materials, with comparison to previous findings 
whenever possible and predictions for future measure- 
ments. The implications of our model and our results are 
discussed in Section V. We end with a summary (Section 
VI) and with ideas for future work. 



II. EXPERIMENTAL MOTIVATION 

Broadening of the temperature-dependent permittivity 
in thin films near the paraelectric-ferroelectric transition 
is reported by several experimental groups 
this observed smearing, accompanied by an overall re- 
duction in its magnitude, is more pronounced with de- 
creasing film thickness. Careful measurements on free- 
standing ferroelectric lamellae yield bulk-type dielectric 
responses, suggesting interfacial effects as the source of 
these thickness-dependent effectsi^ A second related ob- 
servation is that there is a clear separation of temperature 
scales associated with the onset of reversible spontaneous 
polarization and the maximum of the dielectric constant 
in thin ferroelectric filmsi^ 
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Figure 1: Relative permittivity (dielectric constant) as a func- 
tion of temperature (e(T)) for {Ba,Sr)Ti03 (BST) films of 
different thicknesses grown on SrRuOz/MgO (SRO) sub- 
strates^^ where the emax (T) at temperature Tmax is indicated 
with an arrow. 



In Fig. [T] we display relative permittivity measure- 
ments done on {Ba, Sr)Ti03 (BST) thin films grown on 
SrRuOs (SRO).-^ The measurements show suppression 
of the relative permittivity with decreasing film thick- 
ness. As the film thickness decreases to 340 nm, the 
temperature associated with the maximum of the permit- 
tivity (Tmax) appears to saturate at about 300 K, with 
Tmax for the two thinnest films not being clearly discern- 
able. As we will discuss in Section IV, this is consistent 
with the prediction of our model that there should be a 
maximum in Tmax{l) at a thickness I — I*] from the data 
in Figure [T] we estimate I* < 340 nm in BST films grown 
on SRO. 

In ferroelectric films, in contrast with their bulk 
counterparts, there is an observed distinction^ between 
Tmax and Tferro, the temperature where polarization 
becomes switchable. This separation of temperature- 
scales and the permittivity broadening discussed ear- 
lier are both features characteristic of dielectric behav- 
ior in an applied bias fieldlS; the latter could be due 
to a real charge distribution or it could result from an- 
other physical effect^" that breaks the symmetry P 
—P. For example it has been noted- that flexoelectric 
coupling, known to increase near a ferroelectric transi- 
tion,^^ implies a spatially- varying effective field term due 
to the underlying inhomogeneous strain^. The result- 
ing phenomenological model successfully reproduces key 
thickness-dependences of the dielectric properties.— Here 
we ask whether these experimental trends are indeed 
proof of underlying strain inhomogeneities, or whether 
they may be consistent with another strain profile. 




Figure 2: A schematic of the planar ferroelectric capacitor 
under consideration, with the key length-scales and regions 
clearly demarcated. Note that the mismatch defects are seg- 
regated in a buffer layer of thickness Ib and that the polariza- 
tion and strain are homogeneous in the majority of the film. 
Incomplete charge compensation at the ferroelectric-electrode 
interfaces results in a residual depolarization field, as shown. 
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Figure 3: Schematic of the segregated strain model; here the 
elastic defects reside within Ib (<< I) of the film-substrate 
interface so that the majority of the film is homogeneously 
strained. 



III. LANDAU THEORY 

We model the thickness-dependent properties of 
strained ferroelectric films using a Landau approaches' 
where all misfit dislocations are assumed to reside within 
a thin buffer layer of size Ib << I where I is the film 
thickness (e.g. Fig. [21); this is in contrast to other ap- 
proaches where these defects are assumed to be roughly 
uniformly distributed within the film.^ Elastic relaxation 
then occurs so that there is homogeneous strain in the 
film except for the buffer layer (e.g. Fig. [3]). Recent 
X-ray diffraction experiments^S, support previous sugges- 
tions^'' that the in-plane film strain, ui, decreases with in- 
creasing thickness (/) of the overall films; furthermore the 
observed thickness-dependent strain fits an exponential 
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form^'^'^^ very well. Therefore in our segregated strain 
approach, thin films experience homogeneous film strain 
that decreases exponentially with their overall thickness 
I; more specifically we model this thickness-dependent 
in-plane film strain as 



-l/lc 



Um = (b- a)/b, 



(1) 



where Ic is a characteristic length-scale of the strain re- 
laxation, and a and b refer to the in-plane lattice param- 
eters of the film (bulk value) and the substrate, respec- 
tiveljsi^. We emphasize that our values of Ic are deter- 
mined from X-ray diffractio n^'^i^^ for films of a range of 
thicknesses. We note that these values of Ic are compa- 
rable to the film thicknesses of interest {Ic ~ I), result- 
ing in a non-trivial thickness-dependence of the strain- 
related properties. More generally, we remark that the 
segregated strain approach described here has been used 
in modelling epitaxially strained superlattices^^ with re- 
sults that are consistent with experiment ,-22. 

In Fig. [2] we display a schematic of the planar ferroelec- 
tric capacitor under consideration with the length-scales 
involved. More specifically we assume a single-domain 
ferroelectric film with the (uniaxial) polarization in the 
z-direction, normal to the film-substrate interface. Phys- 
ically we expect a build-up of free surface charge at the 
ferroelectric boundaries which, if uncompensated, results 
in a depolarization field Ed- In practice such depolar- 
ization effects are significantly reduced by metal elec- 
trodes that provide charge compensation; however their 
non-ideal nature means that some residual depolarization 
field remains. Because Ed is proportional to y, where 
Ic is the screening length of the electrodes, its impor- 
tance increases with reduced film thicknessii'^'^!^ and 
becomes significant^i^^ for / ~ lOOnm; a term in our 
Landau expansion will be included to account for these 
depolarization effects. 

The defect buffer layer is also displayed in Fig. [2] and, 
as discussed above, we assume homogeneous (but pos- 
sibly relaxed relative to the substrate) strain in the re- 
mainder of the ferroelectric film which is controlled by 
the substrate-film lattice mismatch. Following a pre- 
vious analysis of epitaxially strained films, ^ the stress 
tensors azz = ctxz = fyz = 0, because there are no 
tractions acting on the top film surface. For the spe- 
cial case of a (001) ferroelectric film grown on a cu- 
bic substrate, the strain Uxx = Uyy = ui and u^y — 
since the angle between the two lattice vectors remains 
unchanged {9 = We consider here film tetragonal 
symmetry with Uzz 7^ and u 



yz 



ixz = 0. These 
mixed mechanical boundary conditions associated with 
two-dimensional clamping mean that the standard elastic 
free energy G(P, tr) cannot be used to find the equilibrium 
properties of these systems; instead a Legendre transfor- 
mation, G'{P,ui) = G{P,(T)'\-ui{axx + cryy) + UzzCrzz, to 
a modified thermodynamic potential must be performed 
in order to study the equilibrium properties of the con- 
strained film^ 

We are almost ready to write an expression for 



G'{P,ui) and to calculate observable thermodynamic 
quantities. As discussed earlier, the experiments suggest 
a term in G that breaks the symmetry P — P; this is 
achieved by linear coupling of P to an external electric 
field Eext and/or to an effective bias fieldi^ which we take 
to be of the form 



(2) 



where 1^ ~ Ic- We note that the thickness-dependence 
of Wi is included to model the increased smearing of the 
dielectric susceptibility with decreasing I of ferroelectric 
films.'^" At present we will treat Wi phenomenologically, 
and will defer discussion of its exponential decay and its 
possible origins to Section [V] 

Putting all these elements together, we begin our phe- 
nomenological study with the free-energy expansion 



G(P,CT,T) 



ia(r)p2 + i^P^ - (Wi + Ecxt)P 

— Qll<7zzP^ — Quio'xx + <^yy)P^ 
-7;S1i{(jIx + ^ly + <^lz) - Sl2{(Jxx(Jyy) 



-Sl2Crzz{<^xx + <^yy) — ^■^ii'^lyi 



(3) 



where a{T) = a{T) + a^; a{T) = (3{T - Tbuik), Tbuik is 
the bulk transition temperature, ad is discussed below, 
and f3 and 7 are Landau coefficients; here Qij and Sy are 
the electrostrictive constants and the elastic compliances 
at constant polarization respectively. The depolarization 
field contributes to the free energy through the coefficient 
&{T) in Eq.@?&i^ 



ad 



(4) 



where /g is the screening length of the electrodes, and eg 
and Ee are the electric permittivities of the vacuum and 
the electrodes respectively. 

The mechanical conditions in the film are dG/daxx — 
dG/d(Tyy = —ui, dG/d<Jxy = and dG/dazz = —UzA- 
Solving for the in-plane stresses, one finds that axy = 
and axx = Cyy = cr, where the applied in-plane stress a is 
eliminated by the in-plane strain ui{a). This procedure, 
together with azz = 0, leads to 



G'{Pi,ui,T) = 



Sll 

1 



S12 



7i^*i{T)Pf 



-l*Pt -{Wi+Ecxt)Pu 



(5) 



where 7* = 7+ ^'^''^ anda;(T) = a{T)-ui^^^+ad. 
We note that we explicitly refer to the ^-dependence 
of the polarization {Pi), which here results from the 
thickness-dependence of the strain (u/), the bias field 
{Wi) and depolarization field term {ad), consistent with 
observation.^ One can express the out-of-plane strain 
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(uzz) through its dependence on the out-of-plane polar- 
ization (Pi) and the in-plane strain (it;) as 



2su 



Sll + Sl2 



-Ul- 



-{Qii- 



2s 



12 



Sll + Si2 



!12 



)P?- (6) 



Table I: Landau parameters for BST^ and STCM (in SI 
units). We use the T > lOOK values for a{T), except in 
calculations in Fig. [15] where we interpolate the T > lOOK 
and the T < 50K values to 75K where the two a{T) functions 
cross. 



The Curie film temperature T* refers to the 
paraelectric-ferroelectric transition at zero total field, 
Ef = Wi + Eext = 0. It increases with an applied com- 
pressive strain {Q12 < 0) 



T, 



bulk 



ui- 



4Q 



12 



I. 



Sll + S12 



(7) 



but has a decreasing component for very thin films due 
to depolarization effects. The dielectric susceptibility is 



Xi 



(8) 



which diverges at T* if the spontaneous polarization 
Pi — > 0; that can only occur if Ef' = (see Eq. [10] be- 
low). The dielectric susceptibility is observed to diverge 
for bulk systems^ and for free-standing filmsi^ In gen- 
eral Ef^ ^ for ferroelectric capacitors and P/ at T* , 
so that xi ^^^s a finite maximum at a temperature Tmax 
defined by ^|t=t„„^ = 0. We note that this condition 
combined with the expression for xi above yields 
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7(r) [io«] 


Qii 


Q12 


Sll +S12 [10-^2] 


BST 


4(796 + 2.16(r- 273)) 


0.110 


-0.0430 


5.6 


STO 


1700 


0.066 


-0.0135 


3.0 



fihii 


a(T) [10'^] 


BST 
STO (T > lOOA") 
STO (T < 50K) 


9.1(r- 235.0) 
7.06(T - 35.5) 
263.5 (Coth[42.0/r] - 0.90476) 



where 



Table II: Film parameters: effective field Wo and compressive 
strain Urri^ with associated length-scales Iw and (see 
Eqs. [1] and [2]); also the values for screening length and 
the relative permittivity of electrodes are shown. 



film 


substrate 


Wo [kV/cm] 


Iw [nm] 


u„, [%] 


Ic [nm] 


L [nm] 
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BST 


SRO 


400 


300 


-0.50 


300 


0.023 


1.0 


BST 


PSS 


450 


300 


-0.77 


300 


0.400 


1.0 


STO 


LSAT 


400 


300 


-0.90 


300 


0.023 


1.0 



Pl{T) 



dPi{T) 



dT 



67* 



(9) 



7^ = 



277*3 



47* 



(13) 



where we see that T^^^, differs from T* and depends on 
film thickness via the polarization; this equation gener- 
ally has to be solved numerically to obtain T^^^ once the 
expression for Pi{T) has been determined. 

The condition for finding the system in its equilib- 
rium state is ^ = 0. The spontaneous polarization Pi 
emerges then as the solution(s) to the following cubic 
equation 



and the number of polarization solutions is determined by 
the sign of TZ so that the single nonswitchable Pp corre- 
sponds to 7?. > 0. Therefore the transition temperature 
Tferro between nonswitchable and switchable polariza- 
tion occurs when TZ = Q leading to the expression 



T-. 



ferro 



/? V 4 



1/3 



(Eh 



-T\2/3 



(14) 



anT)Pi+rPi' = El 



(10) 



where for Ef' ^ Wi + Eext + 0, we have to be careful 
to distinguish between the paraelectric (nonswitchable) 
polarization Pp and the ferroelectric (switchable) polar- 
ization P/; here switchable refers to the fact that there 
are multiple solutions for the polarization that can be 
accessed by application of a finite i?ea:t- There are three 
solutions to the equation pO]) 



1/3 



.27 



27* 



1/3 



(11) 



and 



Pf = -P,±i^ 




At this temperature, the paraelectric solution becomes 
an unstable extrcmum. 

In general, the three temperature-scales Tferro, Tmax 
and T* differ as indicated in Fig. [5] We note that for very 
thin films (< 60 nm), there is suppression of all three 
temperatures due to depolarization effects. We also re- 
mark on the presence of a maximum in Tmax{l) that has 
already been alluded to in Section II; this feature will be 
discussed in more detail when we apply this phenomenol- 
ogy to specific materials and substrates. 



IV. RESULTS 

In this section, we calculate dielectric properties for 
two specific materials, (Bao.5SrQ,Q)TiO^ (BST) and 
SrTiOz (STO). Our study of BST films allows us to make 
direct comparison between our calculated properties and 
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the experiments (Fig. [T]) that motivated the inhomoge- 
neous strain scenario.^ In order to explore different pa- 
rameter regimes, we study tliese films on two distinct 
substrates, SrRuOs (SRO) and Pt/Si02/Si (PSS); here 
we note that the latter is a hypothetical case since to date 
epitaxially-grown single-crystal films of BST on PSS have 
not yet been realized. 

We also apply our phenomenological treatment to STO 
films that are known for their coherence; this is achieved 
by highly controlled growth conditions that inhibit de- 
fect formation and thus inhomogeneous strain effects are 
not expected 1^ To our knowledge, there do not exist pub- 
lished high-resolution dielectric measurements of strained 
STO films with polarization normal to the electrode- 
ferroelectric interface. We therefore compare our results 
whenever possible to first-principles calculations,'^^ and 
make experimental predictions for a range of epitaxial 
strain values that could be realized by a variety of sub- 
strates. 

The parameters used in our calculations are presented 
in Tables [T] and [TTl Table [T] indicates the Landau co- 
efficients used for each material. Film-related parame- 
ters, displayed in Table |TT1 are determined from data on 
strain relaxation; the characteristic length 1^ (see Eq. [1]) 
from the lattice constant measurements lies somewhere 
between 200 and 300 nm^^i^l. In order to make compari- 
son with the inhomogeneous strain model scenario, which 
uses a characteristic length scale 300 nm^, we keep this 
Ic value in our calculations. Data on electrode screening 
lengths He) is le of SrRuOs (SRO) 0.23 A^^, and we set le 
of {LaAl03)o.29 X (5^^^0.5^00.503)0.71 (LSAT) to be the 
same value, since LSAT is dominated by S'rAZo.5Tao.503 
that is very similar to SRO. We choose the screening 
length of Pt/Si02/Si (PSS) to be = 4 A, which is 
expected to be larger than l^. in metallic SRO due to the 
presence of semiconducting silicon. We display strain Um 
value for three different substrates, BST on SRO,^ BST 
on PSS {assT = 3.95 A, bpt = 3.92 A) and STO on 
LSAT,^ The bias field is set to be Wq = 4.0 x 10^ V/m 
for BST films on SRO, a value that is comparable to that 
of applied external fields in related BST dielectric mea- 
surements)^ We keep the same Wq for STO films and 
a slightly different one (Wo = 4.5 x 10^ V/m) for BST 
films on PSS. The bias field Wi is treated phenomeno- 
logically (see Eq. [2]) and we emphasize its crucial role in 
modelling key features of the dielectric properties of fer- 
roelectric films as will be discussed in more detail shortly. 



A. BST 

In Fig. m we display the three temperature-scales as 
a function of thickness for BST films on SRO elec- 
trodes with Eext — 0: Tmax, T* and Tf^rro- Because 
EJ = Wi ^ 0, there is a clear separation of the three 
temperatures; for Ef = 0, they collapse onto T* (com- 
pare Eqs.[7l[9land[T4|l. Therefore an estimate of Wq can 
be obtained from the difference of T^ax and T*, where 
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Figure 4: Thickness-dependence of tlie three distinct 
temperature-scales Tf^rro (solid), Tmax (dashed) and T* (dot- 
ted line) with Eext = in the segregated defect model de- 
scribed in the text; here Landau coefficients for BST on SRO 
(see Tables HI and HIl in Section IV) have been used and Tmaxi}) 
is noted to display a peak at I* = 60 nm. 



the latter can be expressed in terms of (experimentally 
accessible) Landau parameters and by the value of the 
misfit strain (Eq. [7]). The magnitude of the obtained 
bias field, Wq — 400 kV/cm, is in rough agreement with 
experiment: the experimental temperature T^]^^ = 250 
K for a 950 nm film in Fig. [1] is close to the calculated 
T^^l = 268 K in Fig. H and T,^^^ = 280 K for the 660 
nm film in Fig.[T]is also reasonably close to the calculated 
T^i''^ = 304 K in Fig. H 

The temperature Tmax displays a peak, Tmax — 519 
K at r = 60 nm, and it decreases with increasing thick- 
ness and asymptotically approaches T*, the bulk tran- 
sition temperature. As previously noted in Section II, 
for I < I* , Tmax decreases with decreasing Z; by contrast 
for I > I*, the trend is consistent with that displayed in 
Fig. [T] More generally, the behavior of Tmax depends on 
that of the polarization as a function of temperature (see 
Eq. [5]) . Both the strain and the bias field make the po- 
larization decrease with increasing film thickness (Z > 60 
nm) as will be discussed shortly, and this monotonic be- 
havior makes Tmax decrease as well. The observed peak 
in Tmax is driven by the depolarization field contribution 
in our model. Tmax becomes suppressed for very thin 
films {I < 60 nm) and approaches zero just as does the 
polarization at these film thicknesses (see Fig. [7]). The 
value of I* = 60 nm is determined by the strength of 
the depolarization contribution, and thus is expected to 
depend on the electrode/substrate material. As already 
noted, we use a longer 1^ for BST films on PSS than for 
BST on SRO. Because the depolarization field contribu- 
tion to our free energy is proportional to y, (see Eq. H]), 
the resulting I* will be larger for BST on PSS than on 
the SRO substrate. For BST on PSS, we calculate that 
the peak in Tmax occurs at I* ~ 300 nm in contrast to 
the value of /* ~ 60 nm for BST on SRO displayed in 
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Figure 5: Calculated relative permittivity as a function of 
temperature for BST films on SRO substrates of thicknesses 
I = lOfim, 950, 660, 340, 280 and 175 nm with Ee^t = 0. 
The highest permittivity corresponds to the thickest film; the 
divergence for 10/im film at the bulk Curie temperature 235 
K is indicated. Reduction in the permittivity for thin films is 
observed; the peak of the permittivity shifts towards higher 
temperatures in agreement with Fig. [1] 
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Figure 6: Calculated relative permittivity as a function of 
temperature for BST films epitaxially grown on PSS sub- 
strates of thicknesses I — 40, 68, 150, 300 and 580 nm with 
Esxt = 0. The highest permittivity corresponds to the thick- 
est film. Reduction in the permittivity for thin films is ob- 
served, and the peak of permittivity shifts towards lower tem- 
peratures. 



Fig.il 

The theoretical temperature T* in Fig. [4] refers to 
the paraelectric-ferroelectric transition at zero total field, 
Ef ^ Wi + Eext (see Eq. [T]). It has a peak at 80 
nm with a maximum value of T* — 329 K and then de- 
creases with increasing film thickness due to the strain 
relaxation term, ui ~ e^'^'"^. It reaches its bulk transi- 
tion temperature value (Tbuik = 235 K) for thicker fully 
relaxed films: T* = 238 K for 1000 nm film. Due to the 
depolarization contribution (see Eq. [3]) , T* is suppressed 
for very thin films and eventually reaches zero. 
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Figure 7: Calculated nonswitchable polarization P[;^] (see 
Fig. [4]) for BST films on SRO substrates of various thick- 
nesses (I) at temperature 300 K and with E^xt = 0; polariza- 
tion as a function of temperature for films of thickness 175 
(solid) and 620 nm (dashed line) are shown in the inset. 



The temperature Tf^rro that separates switchable and 
nonswitchable polarization regimes increases for increas- 
ing / and saturates when it reaches the bulk transition 
temperature. In Eq. [TJ] we see that there is competi- 
tion between thickness-dependent contributions due to ui 
and to Wi in Tferro] therefore whether this temperature 
increases or decreases with increasing film thickness de- 
pends on their relative magnitudes. For example, T/erro 
is shown to increase with increasing I in Figure 21 reflect- 
ing the dominance of the bias field contribution in this 
particular case. We note that the switchable (ferroelec- 
tric) regime is lost as the films become very thin: e.g. 
ferroelectricity vanishes at I ~ 100 nm at fixed tempera- 
ture 60 K. Tferro reachcs zero at a critical film thickness 
{^crit ~ 35 nm here) and films with I < l^rit remain in the 
nonswitchable polarization regime at all temperatures. 

In Fig. [5] we display our calculated relative permittivity 
e(r) (e = 1-f x) for BST films on SRO substrates with 
Eext = 0. Reduction in e(/) as a function of decreas- 
ing thickness is observed for both theory and experiment 
(Fig.[T|); e,nax decreases and T^ax increases as I decreases. 
Favorable comparison of the calculated relative permit- 
tivity e™^^ can be made to its measured analog^^ ^max 



shown in Fig. [TJ e"^]^^ - 
1900 at T, 



with e^^P 



2035 at T,nax = 268 K matches 
^(jj. - — ^',aax ~ 260 K for the 950 nm film, 
and e^i^^ = 1044 at T^ax = 304 K matches with the 
e^P^ « 1100 at Tmax ~ 290 K for the 660 nm film. 

We also predict the thickness-dependent dielectric be- 
havior of e(T) for BST films epitaxially grown on PSS in 
Fig. [HI For films of decreasing thickness where I < I* & 
systematic reduction in e is observed and the peak in the 
permittivity is shifted to lower temperatures. 

Next we present the calculated nonswitchable polar- 
ization in Fig. [7] at temperature 300 K with E^^xt — for 
BST films on SRO. The polarization is roughly propor- 
tional to the bias field Wi and its value increases with 
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Figure 8: Transition temperature Tferro as a function of ap- 
plied external electric field Eext for BST films on SRO sub- 
strates of thicknesses 100 (solid), 175 (dashed) and 620 nm 
(dotted line). The corresponding line divides the region into 
switchable and nonswitchable polarization phases for each 
film. 





0.12 




0.08 


c 


0.04 




N 

M 





"o 

Oh 


-0.04 




-0.08 




-0.12 




-300 -280 -260 -240 -220 
External Field [kV/cm] 



-200 



Figure 9: Calculated hysteresis polarization loops for BST 
films on SRO substrates of thicknesses 100 and 175 nm at 
temperature 290 K. 



increasing misfit strain ui (see Eq. [TU)) . Motivated by 
experiment, we have modelled the bias field and misfit 
strain to decrease with increasing film thickness expo- 
nentially {{Wi,ui} ~ e~'/'"=, Ix — {Iwjc}), and therefore 
the polarization also decreases with thicker films. How- 
ever, for very thin films {I < 50 nm), there is suppres- 
sion of the polarization (P) due to depolarization field 
effects and a peak with a maximum value of polariza- 
tion Pmax = 0.2 [C/to^] at ^ = 50 nm develops. The 
temperature-dependence of the polarization for two dif- 
ferent film thicknesses is displayed in the inset of Fig. [T] 
These P{T) curves shows good qualitative agreement 
with analogous measurements.^ In general, the polariza- 
tion decreases with increasing temperature. 

Tferro, the transition temperature separating the pres- 
ence of switchable and nonswitchable spontaneous polar- 



ization, is plotted in the presence of external electric field 
for BST films on SRO in FigM Tferro{Eext) reaches its 
maximum at Eext — —Wi (where Wi = Woe~'/''=), and 
decreases symmetrically about this value in agreement 
with Eq. 1141 we note that the maximum of T/erro de- 
creases for thicker films as anticipated by the zero-total 
field {Ef) results of T* displayed in Fig.d] The tempera- 
ture Tferro at zcro external field matches the behavior of 
Fig. m Switchable and nonswitchable polarization phases 
are marked. 

Calculated hysteresis loops are displayed in Fig. [5] for 
BST films on SRO substrates at T = 290 K with 1 = 100 
nm and I = 175 nm. According to Fig. [8l at this temper- 
ature for these film thicknesses, the switchable polariza- 
tion develops only for certain values of nonzero external 
electric field: E^xt = {-298, -275} kV/cm for 100 nm 
film and Eext = {-229, -218} kV/cm for 175 nm film. 
The width of the hysteresis loops in Fig. [S] is given by 
the above field ranges; it decreases with increasing film 
thickness and shows good qualitative agreement with ex- 
periment.^^ Hysteresis loops are symmetric around the 
point Eext = —Wi and P = 0, where the bias field Wi 
is the thickness-dependent field offset. This field offset 
(specifically, we refer to the shift of the center of the hys- 
teresis loop along the field axis from zero field position) 
becomes larger for thinner films (see Eq. [J) and is tem- 
perature independent. However, the width of the loops 
shrinks as temperature is increased, as shown in Fig. [S] 
Symmetry in the hysteresis loops around P = yields 
zero offset in the spontaneous polarization and therefore 
no associated charge offset within the thin film. We note 
that significant charge offsets are observed in graded films 
with designed polarization and strain gradients42i 



B. Strained STO 

Bulk strontium titanate SrTiO^ (STO) remains para- 
electric down to the lowest temperatures accessible, but 
strained STO films may be driven into a ferroelectric 
phase even at room temperature.^ To our knowledge, 
detailed thickness-dependent dielectric measurements on 
such STO films have been performed only with ten- 
sile epitaxial strain resulting to in-plane polarization.— 
Here we make predictions for the thickness-dependent di- 
electric properties of STO films with compressive epitax- 
ial strains and polarizations normal to the film-substrate 
interface. 

We begin by making direct comparison between our 
results and those of ab-initio studies, ■^^ displayed in 
Figs. IIOIIII Since ab-initio calculations consider uni- 
formly strained films without strain relaxation and with- 
out an effective bias field, we set Wi = for the purpose 
of comparison here. In Fig. \TU\ we present the spon- 
taneous polarization as a function of misfit strain for 



a I 



100 nm STO film at zero total field, Ef 



0. 



Dots in the figure correspond to first-principles calcula- 
tions;^ where the out-of-plane polarization in the ferro- 
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Figure 10: Spontaneous polarization P[C/m^] as a function of 
film strain (mj) for 100 nm STO film at temperatures T = 200 
(dotted), T = 250 (dashed) and T = 300 K (solid line) and 
zero total field, EJ — 0. Dots correspond to ab-initio values 
of the polarization^. 



Figure 12: Calculated relative permittivity as a function of 
temperature for STO films of thicknesses / = 50, 100, 200, 
300, 400, 500 and 600 nm at misfit compressive strain Um = 
—0.9% with E^^t = 0. The highest permittivity corresponds 
to the thickest film. Reduction in the permittivity for thin 
films is observed. 
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Figure 11: Calculated paraelectric relative permittivity as a 
function of external electric field Ee^^t and film strain it; = 
-0.5, -0.4, -0.3, -0.2, -0.1 and 0.0% for 100 nm STO film 
at temperature T = 200 K. The permittivity data are shown 
in the limit of zero bias field Wi in order to make comparison 
with the ab-initio data^^. 



electric tetragonal phase (u/ < —0.75%) for films with 
zero macroscopic field has been calculated. We choose 
the I = 100 nm STO film where we do not expect depo- 
larization effects to be important [l > l*\ see Fig. [T4)) 
for comparison with the ab-initio data. Good agree- 
ment is achieved at temperatures T ~ 250 K; our cal- 
culated curves follow the behavior of the ab-initio dots. 
At lower temperatures, better agreement is achieved for 
less compressive strain, correctly indicating that fully re- 
laxed STO is paraelectric down to zero temperature. 

Continuing our comparison with the results of ab-initio 
calculations,— we display the paraelectric relative per- 
mittivity as a function of external electric field E^xt and 
film strain ui for a Z = 100 nm STO film at T = 200 
K in Fig. [TT] A nonpolar tetragonal phase develops for 
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T[K] 

Figure 13: Calculated relative permittivity as a function of 
temperature for STO films of thicknesses I — 50, 300 and 600 
nm at misfit strain Um = —0.9% (solid) and Um = —3.0% 
(dotted line) with Eext ~ 0. Higher values of misfit com- 
pressive strain shifts the permittivity curve towards to higher 
temperatures where the larger shifts of Tmax are observed for 
the thinner films. 



strains —0.75% < Ui < -1-0.54% according to the ab- 
initio calculations 1^ ^From Fig. [TOl the best fit for the 
polarization just at the paraelectric-ferroelectric phase 
boundary {ui = -0.75%) is achieved at T = 200 K, 
and therefore we choose this temperature to calculate our 
paraelectric permittivity data. We compare our results 
to the ab-initio calculations in the range of compressive 
strain —0.5% < ui < 0.0%: in both cases, the permittiv- 
ity (e) increases with increasing compressive strain; this 
occurs even more rapidly in the range of external fields 
— 150 < Eext < 150 kV/cm, and its magnitude in the 
range 800 < e < 1500 in Fig. [TT] roughly corresponds to 
the ab-initio values 400 < e < 1800^^ We note that the 
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Figure 14: Calculated nonswitchable polarization P[C /m^\ of 
STO films as a function of film thickness (/) for misfit strains 
= -3.0, -2.5, -2.0, -1.5 and -0.9% at temperature 300 
K and with Eext ~ 0. In the inset, the polarization as a 
function of temperature (T) is plotted for 100 (solid) and 500 
nm (dashed line) films at misfit strain Um = —0.9% and zero 
external field. 
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Figure 15: Transition temperature Tf^rro separating switch- 
able (below) and nonswitchable (above the temperature 
curve) polarization phases as a function of STO film thick- 
ness I and misfit strain Um = —3.0, —2.5 and —2.0 % with 
Eext = 0. Tferro bccomes negative for compressive misfit 
strain values smaller than 1.5% for all film thicknesses, re- 
sulting in a nonswitchable polarization regime. 



observed increase of the paraelectric permittivity with 
strain can be understood from Eq. [8l 

Both Figs . [TOl and [TT] indicate good agreement between 
results of our phenomenological model and those of pre- 
vious first-principles calculations, and this provides us 
with confidence regarding the Landau coefficients and 
more generally with the approach described here applied 
to strained STO films. Next we calculate the dielectric 
properties of strained STO films using a nonzero value 
for the effective bias field, Wq ~ 400 kV/cm, that is 
comparable to that used for BST films in Sec IIV Al 

Our calculated relative permittivity as a function of 



temperature and film thickness at compressive strain 
= -0.9% (STO on LSAT) with Ee^^t = is plot- 
ted in Fig.[T3 The permittivity is suppressed for thinner 
films and its maximum is shifted towards higher tem- 
peratures, displaying similar trends for both e and T^ax 
as for BST on SRO in Sec IIV Al T^ax increases with 
decreasing /, but develops a peak at Z = 60 nm and is 
again suppressed for very thin films due to depolariza- 
tion effects. The screening length of LSAT is comparable 
to If. for SRO, and thus so is /*. The magnitude of e 
is also comparable to that in BST films in Fig. O it re- 
sults from similar values of the Landau coefficients and 
the value of the compressive strain in both films (see 
Table m. To our knowledge, there exists only one pub- 
lished dielectric measurement on strained STO with the 
polarization normal to the film-substrate interface; this 
experiment, performed on a, I — 50 nm STO film grown 
on LSAT (u™ = -0.9%)S. yields e^P^ « 400, that is in a 
good agreement with our calculated maximum value of 
permittivity, e™^^ = 364 for this film/substrate combi- 
nation. We investigate the effect of compressive strain on 
the relative permittivity in Fig. [T31 We observe that in- 
creased compressive strain shifts the permittivity curve 
towards higher temperatures and larger shifts of T^ax 
occur for thinner films. 

We display the nonswitchable polarization as a func- 
tion of film thickness and misfit strain at T = 300 K and 
Eext = in Fig. [131 Again, the polarization shows similar 
behavior as in BST films in Sec. IIV Al it decreases with 
increasing film thickness and is suppressed for very thin 
films (< 50 nm) due to depolarization effects. Increas- 
ing the strain results in higher polarization, in agreement 
with Fig. 1101 this time however the polarization values 
change due to the effect of a nonzero bias field. In the 
inset, the temperature-dependence of the nonswitchable 
polarization is plotted in I = 100 and I = 500 nm STO 
films at misfit strain Um = -0.9% (STO on LSAT sub- 
strate). The polarization decreases with increasing tem- 
perature. 

Unstrained bulk STO remains paralectric down to zero 
temperature. However, as previously noted, a ferroelec- 
tric regime occurs for strained STO films. ^ We plot the 
transition temperature Tferro, separating switchable and 
nonswitchable polarization regimes, as a function of film 
thickness and misfit strain with Eext = in Fig. 1151 
We predict a ferroelectric phase to occur for compres- 
sive strains larger than 2.0%, and note that ferroelectric- 
ity is recovered here for the thinnest STO films, as one 
goes from a thick-film nonswitchable regime to a thin-film 
switchable one (at fixed T); e.g. ferroelectricity emerges 
at —2.0% strained films for thicknesses, 30 < / < 160 nm. 
This is distinct from the behavior previously described 
in BST films, where ferroelectricity is lost by making 
films thinner (see Fig. 3]). Tferro indicates a maximum at 
I* « 60 nm, and this peak is due to depolarization effects 
(see Eq.U]) that reduce the transition temperature to zero 
for the thinnest films. Tferro decreases in thicker films 
{I > 60 nm) for values of compressive strain Um > 2.0% 
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Inhomogeneous Strain Model 
Strain u(z) ; 




Figure 16: Schematic of the inhomogeneous strain model^- 
where the effective bias field W{z) is spatially-dependent due 
to flexoelectric coupling; here the characteristic length-scales 
associated with the strain relaxation (Ic) and the buffer layer 
(Ib), where the elastic defects reside, are assumed comparable 
with the overall film thickness (l). 



as shown in Fig. 1151 However, it increases with increas- 
ing film thickness for small values of compressive strain 
{um ^ 1.5%), similarly to EST films on SRO substrates 
(with Um = —0.5%) in Fig. [4] We note here that although 
Tferro incrcascs, it has negative nonphysical value for 
these low strain values, and thin films remain in the non- 
switchable regime down to zero temperature. As noted 
previously, the thickness-dependence of this temperature 
scale arises from competition between strain and bias 
field contributions (see Eq. [T3]); the former dominates 
for large enough mistmatch strains (um ^ 2.0%), and in 
this case Tferro decreases with increasing film thickness. 

V. DISCUSSION 

Next we explore the implications of our results and 
the origin of our model assumptions. We begin with a 
general discussion of the effective bias field in epitax- 
ial perovskite oxide films. Both the inhomogeneous (see 
Fig. [in]) and the segregated (see Fig. [3]) strain mod- 
els describe the thickness-dependent dielectric proper- 
ties of ferroelectric films consistently with experiment, 
and thus further measurements are required to determine 
the presence/ absence of underlying inhomogeneous strain 
throughout the film. These models each have effective 
bias fields, one that is spatially varying^ and the other 
that is uniform, and it is exactly this feature that we ex- 
ploit in a proposed benchtop experiment to distinguish 
these two scenarios. 

An effective bias field breaks up-down symmetry at all 
temperatures. In a film above the zero- (external) field 
{Eext = 0) ferroelectric transition temperature Tferro, 
this results in a nonzero macroscopic polarization in zero 



external electric field (see Eq. [T0|) ; for the sake of com- 
pleteness, we note that the Curie temperature T* refers 
to the paraelectric-ferroelectric transition at zero total 
field {Ef = Wi + Eext — 0). While this polarization can 
vary with temperature, making the film pyroelectric, it 
should not be confused with a ferroelectric spontaneous 
polarization. The correct distinction between polar and 
ferroelectric films is made on the basis of switchability, 
as determined, for example, through a hysteresis mea- 
surement. In a nonswitchable polar film, there is only 
one locally stable polarization state with Eext = 0, and 
the system will show dielectric behavior in a finite ex- 
ternal electric field. On the other hand, if there are two 
(or more) stable polarization states for the system with 
Eext = that can be switched by application of an ex- 
ternal electric field, then the film is ferroelectric, and it 
will exhibit a characteric ferroelectric hysteresis loop. An 
effective bias field Wi will lead to a lateral shift in the fer- 
roelectric hysteresis loop, which can be used to determine 
the associated effective bias voltage across the film. Field 
offsets have been indeed experimentally observed in 100 
nm PbZro.2Tio.sO3 (PZT) thin films.^? We also remark 
that an effective bias field will make the two polarization 
states symmetry-inequivalent with different energies. 

The transition from nonswitchable to switchable fer- 
roelectric phases usually occurs as a function of de- 
creasing temperature at Tferro- The detailed thickness- 
dependence of Tferro depends on material-specific pa- 
rameters, as can be seen in Figs. |3] and [12] for EST 
and STO thin films. More generally the Z-dependence of 
the dielectric properties enters via the strain relaxation 
(Eq.[T]), the bias (Eq.[2]) and the depolarization fields. For 
the strain relaxation, an exponential decay on a charac- 
teristic length scale of several hundred nanometers was 
observed experimentall y ^^1^^ In our model, we assume 
the same exponential decay for the magnitude of the uni- 
form effective bias field. These two quantities determine 
the thickness-dependence of the quantities of interest in 
all but the very thinnest films, where the depolarization 
field term dominates, strongly suppressing Tferro, the po- 
larization, and Tmax- In the case of the temperature 
Tferro (sce Figs. HI and [T^ ■ the strain and effective bias 
contributions opppose each other (Eq. I14p and depending 
on their relative strengths, 7/erro(0 increases (EST case 
where Wi dominates) or decreases (situation for strained 
STO where ui is greater) with increasing film thickness 
I. 

A direct consequence of the strain contribution to 
Tferro IS that wc predict that ferroelectricity can be 
strengthened as the films get thinner (ui increases with 
decreasing I), resulting in a transition from a nonswitch- 
able polar phase to a ferroelectric state below a critical 
thickness Ict', more specifically, for STO measured at 
100 K and compressive strain —2.5% (see Fig. [TJ]), the 
critical thickness below which ferroelectricity appears is 
IcT — 200 nm. This runs counter to the usual notion 
that ferroelectricity is suppressed as the film thickness de- 
creases, disappearing below a critical thickness; we note 
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that would be the case for BST on SRO (here Wi dom- 
inates the expression for Tff,rro) where our results dis- 
played in Fig. |4] indicate a critical thickness of about 100 
nm at temperature 60 K. This reentrant ferroelectricity 
as a function of decreasing / should be readily observable 
in an appropriate experiment for strained STO films with 
the polarization normal to the film-substrate interface. 

Both the strain and the effective bias contributions act 
to decrease TmaxQ-) and Pi (related by Eq. [9]) as a func- 
tion of increasing film thickness {I) (see Figs.l4l[7landfT4|). 
By contrast, the depolarization contribution suppresses 
Trnax{l) SLud Pi with decreasing I. The dominance of this 
depolarization term explains the observed shift of the 
peak of Tmax{l) to higher values of I* for semiconduct- 
ing substrates (e.g. PSS) (see Fig ^ that have longer 
screening lengths (le) than their metallic counterparts 
(e.g. SRO); here we recall that the depolarization contri- 
bution to the free energy expansion is ad ~ le/l (Eq. [4]). 
We note that the thickness-dependent effect of the de- 
polarization field on the relative permittivity has been 
noted before-^*' with a similar term, a^^ = a/l, where 
boundary conditions for the spatially-varied polarization 
are proposed that incorporate the effects of a symmetry- 
breaking effective field. In this previous approach^ a 
then is a boundary-related characteristic length. Since 
the thickness-dependence in both treatments is the same, 
one obtains similar results for the relative permittivity 
with appropriate choice of these length-scales {le and a) 
although their physical origins are different. Here we 
have extended this treatment to address the thickness- 
dependence of other dielectric properties (e.g. polariza- 
tion) as well, and we note that the previous inhomoge- 
ncous strain approach^ did not include such depolariza- 
tion effects for thin ferroelectric films. 

The smearing of the sharp peak in the temperature- 
dependent dielectric response (Figs. [5l [6] and [12]) in zero 
electrical field {Eext = 0) is a signature of the presence of 
a finite effective bias field W; this point has been much 
discussed previously both here (see Eqs. [51 [TUj) and by 
others . Wi assumes larger values for thinner films 
(see Eq. [2]) and so pushes the permittivity to smaller val- 
ues in thinner films (see Eq. [S|), in accordance with ex- 
periment (Fig. [T]). In the bulk limit, both the strain and 
the bias field vanish and bulk behavior of the dielectric 
properties is recovered. 

While there is general agreement that effective bias 
fields play an important role in the properties of per- 
ovskite thin films their specific origins and their 
spatial natures in the films are less well understood. For 
example, an effective bias field can be produced by a spa- 
tially (z)-dependent strain via a flexoelectric effectj^ we 
will refer to this as the inhomogeneous strain model and 
it is schematically depicted in Fig. 1161 In this scenario, 
misfit dislocations are distributed roughly isotropically 
throughout the film and produce strain gradients. By 
contrast, in the segregated strain model presented here, 
the elastic defects are concentrated in a thin buffer layer 
near the film-substrate interface (see Fig. [2] and Fig. [3]) . 



Benchtop Probe for the Segregated Strain Model 




Figure 17: Schematic of a benchtop probe to test for the seg- 
regated strain model: the field offset (Wi) in the observed 
hysteresis loop can be tuned to zero by application of an 
electric field E^xtil) = ~Wi; in this case, the relative per- 
mittivity sharpens since the net (thickness-dependent) field 
ET = E:,t{l) + Wi=0. 



This buffer layer itself breaks the up-down symmetry of 
the film, which then results in a nonzero effective field. 
To be more specific, it may be that an edge dislocation 
in the buffer layer produces a local polar distortion. This 
makes the buffer layer polar and produces a field in the 
uniform film. We then expect the strength of the effective 
bias field to be related to the areal density of dislocations, 
and thus to the magnitude of the homogeneous strain in 
the film. 

As we have shown in Sec. IIV Al using the segregated 
strain model (see Figs [2| and [3]) we recover the thickness- 
dependent dielectric properties of BST films consistent 
with those measured^'^ and calculated using flexoelec- 
tric effectsS. within an inhomogeneous strain scenario 
(see Fig [TB)) . Therefore, in order to determine which of 
these two models is applicable to a particular film, addi- 
tional experimental characterization is necessary. The 
presence of strain gradients in the BST films studied 
using flexoelectricity^ was verified by x-ray analyses^, 
though similar studies on different BST films (on dif- 
ferent substrates) indicate the absence of such inhomo- 
geneous strain.— Therefore it is clear that observations 
of thickness-dependent permittivities and polarizations 
are not enough to indicate the underlying strain profile 
of the ferroelectric film. Naturally the presence/ absence 
of strain gradients can be addressed directly by x-ray 
diffraction but this probe may not be always easily ac- 
cessible or practical. 

Here we emphasize that the effective bias fields associ- 
ated with these two strain models are spatially different 
(see Figsll6landl51). and we will use this distinction to pro- 
pose a benchtop experiment to distinguish between these 
two scenarios. In the inhomogeneous strain modelj^ the 
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effective bias field is spatially nonuniform; more specif- 
ically is proportional to the strain gradient via a flexo- 
electric coupling, and thus has an exponentially decaying 
spatial profile within the film (see Fig fTB)) . By contrast, 
in the segregated strain model (see Fig [3]), the effective 
bias field is uniform from the edge of the buffer layer to 
the surface. Assuming that this effective bias field Wi is 
uniform in the majority of the film, one can tune an ex- 
ternal electric field to the right value, E*^f{l) = — Wi, to 
create a net zero-field condition {Et{1) = E*^^{1) + Wi) 
where the specific value of the necessary external field 
would be thickness-dependent. Therefore there would ex- 
ist an external field value E*^f.{l) when the hysteresis loop 
would no longer have a field offset; at this value of the 
applied external field, a sharp peak in the temperature- 
dependence of the relative permittivity should be ob- 
served (Fig[T7)). We emphasize that this must be a unipo- 
lar hysteresis experiment, with single top and bottom 
electrodes; many measurements are performed across two 
top electrodes, namely two series capacitors, where one 
would always be uncompensated. Another probe of the 
spatial uniformity of W would be to measure Tjerro as 
a function of Eg^t for a uniform effective bias field there 
would be a sharp peak in T/erro, as shown in Fig. [8] This 
pronounced peak would not be present for a spatially- 
varying effective bias field W{z) since the latter would 
have varying magnitude in the sample and no particular 
value of applied (uniform) E^xt could completely com- 
pensate for it everywhere in the film; physically we note 
that W{z) could arise from coupling of the polarization 
to gradients in the strain;^ in chemical composition, and 
in temperaturcjiSiiSiii^ There do exist graded ferroelec- 
tric structures where such spatially varying quantities 
are explicitly present by design; here a charge offset is 
often observed in the charge-voltage hysteresis loopSf^ 
and this could serve as an indicator of underlying gra- 
dients in ferroelectric films if a suitable "charge origin" 
could be chosen as a reference. 

We emphasize that we expect different films, with vary- 
ing compositions, substrates and growth conditions, to 
have diverse strain and effective bias field profiles. The 
inhomogeneous scenario may describe some while oth- 
ers may be better modelled by the segregated strain ap- 
proach; still others may exhibit intermediate behavior. 
By carefully monitoring growth conditions, it should be 
possible to control the density and spatial distribution 
of strain-relieving defects; in some case, it may even be 
feasible to kinetically inhibit them to obtain uniform 
coherently-strained films. ^ We emphasize that in each 
case, the strain and effective bias field distributions must 
be carefully characterized for a full interpretation of the 
measured thickness-dependent dielectric behavior, and 
we have presented simple proposals for benchtop mea- 
surements to ascertain the importance of strain gradients 
in the films. 



VI. SUMMARY 

In conclusion, we have demonstrated that a segregated 
strain model (cf. Fig [J) describes the observed thickness- 
dependent dielectric properties of ferroelectric films as 
well as does a previous model of inhomogeneous strain?^ 
Therefore such thickness-dependent behavior is not sig- 
natory of underlying strain gradients, and more measure- 
ments must be performed to determine the strain profile 
in the film. If the effective bias field is spatially uniform, 
it can be compensated by the application of an applied 
external electric field E*^^{1) — —Wi whose value will 
be dependent on the overall film thickness /. Benchtop 
experiments performed with E^xt = ^exti^) '^i^l yield 
bulk-like sharp dielectric responses. However such com- 
pensation will not be possible if the effective bias field 
is spatially varying, since then its effects cannot be can- 
celled by the application of an external uniform field. 

We have compared our results with experiment (BST 
on SRO) whenever possible and have also made predic- 
tions for measurements on strained STO films with out- 
of-plane polarization. Agreement with existing ab-initio 
calculations, when appropriate, has been good. The pos- 
sibility of reentrant ferreoelectricity in strained STO films 
has also been discussed and we hope that this will be ex- 
plored experimentally in the near future. 

Our phenomenological study of planar films suggests 
that their thickness-dependent dielectric properties are 
not indicative of underlying inhomogeneous strain, and 
are consistent with other strain profiles. We view this 
project as the beginning of a broader study of the phys- 
ical consequences of boundary-induced effects in ferro- 
electrics of increasingly complex host topologies. A next 
step is to explore cases where the strain gradients will be 
induced by geometry: examples include curved films and 
cylindrical shells. Because of the coupling between the 
elastic and the electrical degrees of freedom in these sys- 
tems, we expect tunable strain gradients to stabilize novel 
polarization configurations with rich phase behavior, and 
here flexoelectric effects should definitely be investigated. 
More complex host geometries and boundary conditions 
are expected to favor more novel orderings and dielectric 
properties; the possibility of identifying and characteriz- 
ing these features in three-dimensional ferroelectrics on 
the nanoscale could also be useful in the design of future 
ferroelectric memories fi4 
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